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On The Indefinite Snm In Fractional Calcnlus 


James Nixon 


Abstract. We present a theorem on taking the repeated indefinite summa¬ 
tion of a holomorphic function (j){z) in a vertical strip of C satisfying expo¬ 
nential bounds as the imaginary part grows. We arrive at this result using 
transforms from fractional calculus. This affords us the ability to indefinitely 
sum more complicated functions than previously possible; such as holomorphic 
functions of order n E N that have decay at plus or minus imaginary infinity. 
We then further investigate the indefinite summation operator by restricting 
ourselves to a space of functions of exponential type. We arrive at a second 
representation for the indefinite summation operator, equivalent to the first 
presented, and show we have defined a unique operator on this space. We de¬ 
velop a convolution using the indefinite sum that is commutative, associative, 
and distributive over addition. We arrive at a formula for the complex iterates 
of the indefinite sum (the differsum), using this convolution, that resembles 
the Riemann-Liouville differintegral. We close with a generalization of the 
Gamma function. 


1. Introduction 

Indefinite summation is an operation in mathematics that has a long history. 
Many mathematicians have investigated the operator, including Ramanujan [1] and 
Euler [4]. The reasoning behind defining indefinite summation is to give mathe¬ 
matical sense to the analytic continuation of sums of a holomorphic function. As 
in, suppose S'(n) = fiJ) u G N, what happens if we replace n with z € <C 
and require S be holomorphic? There exists techniques for generating the indefi¬ 
nite sum of holomorphic functions and each has its advantages and its flaws, where 
the majority of the results are accumulated using Bernoulli polynomials and the 
Laplace transform. In this paper we will distance ourselves from these methods 
and approach the problem with the differintegral. We will try to explain our best 
the advantages of our technique. 
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The method we introduce for finding the indefinite sum of a holomorphic func¬ 
tion employs a triple integral transform that comes from fractional calculus. Al¬ 
though clunky at first sight, many of the results work out cleanly and with fair 
ease, despite the mess of symbols. The advantage of our method, in comparison 
to Ramanujan’s method, is that we achieve a more localized form of the operator 
on holomorphic functions; whereas Ramanujan’s method requires the function be 
analytic in a half plane of C, we only require it be analytic in a vertical strip of 
C. Also, unlike Euler-Maclaurin and Ramanujan summation, our method of per¬ 
forming indefinite summation only requires knowledge of the function on a vertical 
line in C (or its values at the positive integers), and none of its derivatives are 
used in the formula. Our method also allows us to perform an unlimited number of 
indefinite summations in a single concise formula, which is much more cumbersome 
using Euler-Maclaurin summation or Ramanujan summation. Eurther, our method 

only requires an exponential bound as the imaginary part of the argument grows, 

2 

this implies we can take the indefinite sum of functions of higher order like 
since it has decay at plus or minus imaginary infinity, which is not possible with 
Ramanujan and Euler-Maclaurin summation. 

A large problem that needs to be addressed when investigating the indefinite 
sum is the uniqueness of a particular indefinite summation operator. By example, 
we can always add a one-periodic function to the indefinite sum and still get another 
indefinite sum [1] . For this reason an additional criterion is added to the indefinite 
sum so that it satisfies uniqueness. We do not encounter this problem when defining 
indefinite summation with our transforms because adding a one-periodic function 
to our indefinite sum causes our integral expression to blow up at plus or minus 
imaginary infinity and it no longer converges. This ensures that our indefinite sum 
is unique, although it satisfies moderately stricter restrictions. 

We elaborate further on a series of results involving the indefinite summation 
operator that we can arrive at. We define a space we can perform indefinite sum¬ 
mation in, and from this we develop a second representation of the indefinite sum 
equivalent to the first. This representation only requires knowledge of the indef¬ 
inite sum at the positive integers. From here we are given open ground to talk 
freely about interesting formulae which arrive due to the recursive behaviour of the 
indefinite sum. 

We define a convolution involving the indefinite sum. It allows us to express the 
“differsum” in a formula that resembles the Riemann-Liouville differintegral very 
much. We arrive at these results by only observing the behaviour of the indefinite 
sum at the non-negative integers by using a convenient and powerful lemma. We 
close by first iterating the indefinite sum to an arbitrary complex number in the 
right half plane, and then the indefinite product of the Gamma function. This 
produces a generalized function that interpolates z and r(z) while maintaining a 
recursive pattern. 

With these formalities we introduce the indefinite summation operator. 
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Definition 1.1. If ^ is holomorphic on open 17 C C then an indefinite sum ^ 
of (j) is denoted ^ i® holomorphic on fl and for z, 0 + 1 G 17 satisfies, 

^ + 0(z + 1) = ^ 0 

2 2+1 

We note the strongest result on indefinite summation in complex analysis that 
we know of, Ramanujan’s method, which produces an indefinite sum for any holo¬ 
morphic function (j) of exponential type a < 2^ (as in \4>{z)\ < [1]. So far 

we have not defined the indefinite sum quite rigorously, as issues of uniqueness need 
to be addressed, however, we are satisfied with this formal definition until we bring 
forth our analytic integral expression for the operator. 

We introduce the backwards difference operator. It appears in combinatorics 
and recursion quite frequently. Despite its simple definition, considering it for 
analytic functions proves trickier in comparison to the derivative operator (of which 
it is a discrete analogue). 

Definition 1.2. If ^(z) is holomorphic on open 17 and there is an open subset 
17' C 17 such that zG 17' => z,z—1 g 17 then the backwards difference of 4> on 
17' is given by, 

= Cj){z) - (l){z - 1) 

The backwards difference operator has appeared in many places in mathemat¬ 
ics. The operator has undergone much investigation, but now the concept-like the 
indefinite sum-in the field of calculus and analysis, exists in a niche. The intertwin¬ 
ing property of the indefinite sum and the backwards difference operator is simple 
^V = VX] = for^ = Thus, these operators are inverses of each other. 

Before progressing we introduce a function that is invaluable to our research, 
Euler’s Gamma function. 

Definition 1.3. The Gamma function T is a meromorphic function on the 
entire complex plane with poles at the nonpositive integers, for z G C, 5R(z) > 0 
represented by, 

pCO 

(1.1) r(z)=/ e-H^-Ut 

Jo 

The Gamma function satisfies the functional relationship, zT{z) = r(z -I- 1) 
and converges to the factorial for natural values, r(n-|-1) = n!. We have imaginary 
asymptotics of the Gamma function, which are, 

( 1 . 2 ) \T{a + iy )\^ ?/—>■ ±00 

There exists a plethora of material on the Gamma function and a quick search 
will give papers unlimited. 

Before introducing our results, we present our differintegral, which will be the 
hidden tool we will be applying throughout. For our purposes we only define the 
differintegral centered at zero. 
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Definition 1.4. Let f{x) be holomorphic in the open sector H = {x S C; tt < 
r] < arg(—x) < rj' < tt; 77 , 77 ' G R, 77 ' — 77 tt}. Assume there exists b G 
dependent on / such that for a G 0 < tr < 6 and 0 G R, 77 < 0 < 77 ' we 
have /o“ |/(-e'® x)|x'^ ^ dt < 00 . Then, the Weyl differintegral centered at zero 

of / is represented for z, 0 < 5R(2;) < b, 

x—0 



( 1 . 3 ) 


dx 


piSz 

f{x) = 

x =0 i (z) 


/(-e'®x)a 


dx 


The differintegral as we’ve defined it is holomorphic in z. We can retrieve / 
from its differintegral centered at zero by using Mellin’s inversion theorem. 


1 ncr+ioo 7-5 

^T^lJa-ioo ^dx ^ z:=0 / 

This result shall give us a lot to work with in our arsenal in order to handle 
the problem of producing an indefinite sum. The last operator we introduce is 
the Twisted Mellin Transform. This operator presents a strong relationship be¬ 
tween differentiating a function and performing the backwards difference of the 
transformed function. The operator is introduced in , however we found it inde¬ 
pendently and write it in a slightly generalized form. 


Definition 1 . 5 . Suppose / is holomorphic in the open sector 17 = {x G 
C; —7r/2 < —77 < arg(x) < 77 < 7r/2 ; 77 G R+}. Assume for cr G R'*' in the strip b G 
R"*", 0 < cr < 5 and 0 G R, —77 < 0 < 77 we have dx < 00. 

Then, the Twisted Mellin transform of / is represented for 0 < < b, 


yzf = 


d-^ 

dx~^ 


e"/(-x) 

:0 


The fundamental property of the Twisted Mellin transform is, 

( 1 . 4 ) y.^ = vzyzf 

which follows by. 


dx 




[e"/'(-x)] = 


dx ^ 
d-^ 


dx 


e"/(-x) - 

rc=0 L dx 

a7=0 dx ^ 


X—0 


7 (-x) 


With these introductions we see we have enough to state our analytic contin¬ 
uation of the summation operator. We hope the brevity of our elaboration on the 
results outside of this paper we shall use does not affect its comprehension. 


2. The Analytic Continuation of the Summation Operator 

We introduce a few lemmas on the absolute convergence of certain modified 
Mellin transforms. These lemmas shall follow quickly and are applied one after 
the other to the triple integral transform that will be our method of analytically 
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continuing the summation operator. The first lemma is inspired by the Paley- 
Wiener theorem on the Fourier transform, where they show holomorphic functions 
bounded on a horizontal strip in C by have a Fourier transform bounded by an 
exponential in the same horizontal strip . Our result appears much more fitting 
to our differintegral however. We prove the result using the differintegral, but all it 
takes is some clever re-maneuvering to make it apply on the Mellin transform-which 
we leave to the interested reader. 

Lemma 2.1. Let (j) be holomorphic in the strip 0 < 5i(z) = a < b for b, tr, S 
and let \ lim 3 j(j,) ^o(l){z)\ < oo. Suppose that \'p{z)\ < for a S K, 0 < 

a < 7r/2 and C G K’*'. Then f{x) = ^ is holomorphic 

in the sector a — 7r/2 < arg(a;) < 7r/2 — a and \f{e’'^x)\x'^~^ dx < oo for 
0 G R, a — 7r/2 < 9 < -n— a. 

Conversely, let f be holomorphic in the sector a G U, a — 7r/2 < arg(x) < 
TT12 — a. Suppose for a G R'*" in the strip b G R"*", 0 < a < b and 0 G R, a — 7 r /2 < 
9 < -K12 — a we have \f(e'^^x)\x'^~^ < oo. Then (j){z) = f{x)x^~^ dx is 

holomorphic in the strip 0 < 5R(z) < b and \4>iz)\ < for M G R'*'. 

Proof. We prove the first direction. Since \4>{z)\ < for a — 7r/2 < 

9 < tt/ 2- a and X G R+, /(e*®x) = ^ d^ f is holomor¬ 
phic because 5 ^ d? ^ 2 ¥I la-i^ uniformly in a; 

in the open sector. We also know that |/(e*®a:)| < la-i^ I dz = 

CcrX~'^ where the integral is absolutely convergent again because of the bounds 
on (j) and the asymptotics of the Gamma function (HU. Now we can shift cr in 
the open strip and this does not affect convergence so that |/(e*®a;)| < Cx~'^ as 
a; —>■ 0 for every e > 0 and |/(e*®a;)| < Cx~'^ as a; —>■ 00 for C G R"^ and C = 
sup<,g[,,h]{^ |r(0</'(C)e“*®^| d^}- This ensures that |/(e*®a;)|x'^-i dx < 

00 for 0 < cr < d and a — 'k/2 < 9 < -k 12 — a. 

For the other direction we have 

poo poo 

T{z)(j>e{z) = / /(e*®x)a;^“^ da; = / f{x)x^~^ dx = e~^^^T{z)(j)o{z) 

Jo Jo 

which follows by contour integration. Now take the absolute value to see that 
|0o(-2)| < Co- '|r(z)| where a — 7 r /2 < 9 < ■n/2 — a and Co- = /g°° |/(e®®x)|x‘^“^ dx. 
By the asymptotics of the Gamma function (11.211 . |0o(-z)| < which 

by taking \9\ = 7 r /2 — cr we get \(j)o{z)\ < where (j)o = J) and M G R"*". 

We note we have dropped the term in the inverse Gamma function’s 

asymptotics because it decays to zero as y —>■ ± 00 . □ 

What we have just shown is subtle but very useful. We have not seen it 
proved before and so we decided to show a proof-though we have our suspicions 
a proof is already out there. It shows we have an isomorphism between spaces of 
holomorphic functions (the space / is in to the space (j) is in) which is something we 
can manipulate freely. Our next lemma comes from asymptotics for the incomplete 
Gamma function and is the core of the convergence of our triple integral. 
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Lemma 2.2. Let f{x) he holomorphic in the sector $2 = {x S C; a — tt/2 < 
arg{x) < 7r/2 — a}. Suppose for a S R'*' in the strip b € K+, 0 < a < b and 
0 G R, Of —7r/2 <9 < 7r/2 —a we have dx < oo. Then for \ G R+ 

and 0 < a < b, 

J e^^\f{e^^x)\dx^dt<oo 

Proof. Let us first note that /(O) < C so that J*e^f{x)dx < oo. Perform 
integration by parts: 


pOO nt poo poo 

/ e^“|/(e*^x)|(ixfit= / e^“|/(e*®x)| / e-^H^-^tdx 

Jo Jo Jo Jx 

From looking at we have an asymptotic expansion of the incomplete Gamma 
function in the integral, e~^*t^~^ dt ^ {Xx-a)>‘+^ which implies 

we can find a C G R"*" such that dt < Ce~^^x°'~^. Where this follows 

because we can factor 1/x from the sum and the rest of the series is 0(1) as x —>■ oo. 
This shows, 


poo poo poo 

/ e^“|/(e*®x)| / e-^H'^-Utdx <C |/(e*®x)|x'"-^ dx 

Jo Jx Jo 

This gives the result. 


□ 


We can now give our theorem on the analytic continuation of the summation 
operator. We may not have given enough motivation as to why the integral trans¬ 
forms we present are the ones we are looking for, but the recursive behaviour is 
satisfied and from a familiarity with fractional calculus it seems less questionable 
and more apparent. The expression may seem to have been plucked from the air 
since we have tried to separate the direct attachment to the differintegral, however 
someone who is familiar with the transforms will clearly see that we are conjugating 
(in terms of abstract algebra) the integral operator by the Twisted Mellin trans¬ 
form, where the inverse Twisted Mellin transform is calculated using the inverse 
Mellin transform. 

Theorem 2.1. Let 4>{z) be holomorphic in the strip b G R^, 0 < 5R(z) < b, 
b > 1 and let \ limj}( 2 )_j.o 4>{z)\ < oo. Let \(l){z)\ < for a G R, 0 ^ a < tx 12 

and C G R"*". Then the indefinite sum of (j), f is represented for 0 < a < b and 
0 < 5R(2;) < b by, 

Ri) = Is(/™ 

or written more compactly, 

( 2 . 2 ) Y.f = y,ry-^f 

. do 
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Proof. We first show convergence of the triple integral. This will not be 
too difficult as the bulk of the work was placed in the lemmas before. If f{x) = 
^ la-ix it is holomorphic in the sector a—7^12 < arg(a;) < 7r/2—a 

and satisfies the conditions of Lemma l2.ll so we know, \f{e^^x)\x’^~^ dx < oo for 

a G in the strip 0 < u < &. Therefore / satisfies the conditions of Lemma 12.21 
which implies J* e^\f{e'^^x)\dxdt < oo. Therefore the triple integral 

transform is convergent. Also, g-cos(e)t^(T-i| dx\ dt < oo because 

all of these transforms converge in sectors, so in the last step we are justified in 
applying the Twisted Mellin transform. 

Now we show that this is in fact an indefinite sum of (j). Let us assume that 
1 < 5R(0) < b so that Vz J2z ‘t’ exists and we get, by applying Equation (11.4L 


= Vz 


poo nt 

-tiZ — 1 I X , 


uTT / ® / e."'f{x)dxdt 


1 r d 

r( 7 :) Jo ^ Jo 


^ f dx'j dt 


1 

/ e-H^-^e*f{t)dt 
Jo 


r(z) 


1 


r(^) JO 
= J^iz) 


pOO 

/ f{tr-^dt 

Jo 


Therefore we know that this operator satisfies vS = 1- Let us now show 
it happens on the right as well. We first note that ^ ~ “ 

l))a;“« di = f{x) + f'{x) and that f* e^(f(x) + f'(x)) dx = e*f{t) - /(O), 




1 /■“ /■* 

/ e^{f{x) + f'{x))dxdt 
^ \Z) Jo Jo 

^ \Z) Jo 
J^iz) - /(O) 

(j){z) + A 


Therefore we have shown \/J24‘ = 4‘ S V<?^ = (j) + A (where A G C, and 
this constant is much like a constant of integration, as in g'{t) dt = g{t) + A) so 
that it is natural to say = Sj which is a defining property of ^ and assures 
that YJz + 1) = Sz+i This gives the result. □ 


More advantageously we can show that is exponentially bounded as its 
imaginary argument grows so that we may consider (j). Going on in such a 
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manner we can produce all the natural iterates of the indefinite summation opera¬ 
tor. 

Lemma 2.3. Let 4){z) be holomorphic in the strip b G M’*', 0 < < b for 

b > 1 and let \ limjf(j,)_>o (j){z)\ < oo. Assume \(j>{z)\ < for a S R, 0 ^ a < 

7r/2 and C e R+. Then, | Ez ^1 < for M e K+. 

Proof. Observe the representation of the indefinite sum in m and observe 
that the hnal integral transform applied is the Mellin transform. Observe further, if 
/(^) ~ ^ la-ix df we know / is holomorphic in the sector a — 'Kjl < 

arg(a:) < 7r/2 — a. If a — 71/2 < 0 < 7r/2 — a then, t^~^e~^ f* e^f{x) dxdt = 

/o°° t^~^e~^'^*' Jp dxdt. We have done a substitution through both 

integrals which is justified since all these transforms converge in sectors. We know 
that 9 is arbitrary and the final function does not depend on 9 by considering 
contour integrals. We leave the details to be filled by the reader. These integrals 
are absolutely convergent, again by Lemma \2A \ and Lemma [2?^ By Lemma [2^ we 
also see that 


POO nt 

I / ‘ / e® ^ f (e^^x) dx dt\ < 

do do 


< 


for some K G This shows that. 


— — t cos{9) 


^|/(e*®a;)| dxdt 


pOO 

K \f{e^x)\x’^-^ dx = Ca 
^0 






r(z) 


pOO Pt 

/ / e^"’^f{e^^x)dxdt 

Jo Jo 


so we have | Ez ^1 < ®)l&(^)l which gives us | Ez'/’I for an 

M G R+. This shows the result. □ 


With this we can state our corollary on the fact we can take any natural iterate 
of the indefinite sum and it will remain holomorphic in the same strip. 


Corollary 2.1. If (j){z) is holomorphic in the strip b G 0 < < b 

for b > 1 and \<l){z)\ < for a G R, 0 < a < 7r/2 and C G R+ and 

I limjf( 2 )_,.o (/>(z)| < oo then </> can be indefinitely summed an arbitrary amount of 
times 


(2.3) 


E 


_ pX pX 

...{ntimes)...2^4> = yz / ...{ntimes)... / 
do do 


Proof. Observe Theorem O and Lemma O Theorem o ensures we can 
indefinitely sum (f and then Lemma [231 ensures that that indefinite sum satisfies the 
original conditions of Theorem 12.11 The representation (12.31) follows by induction 
and by applying Lemma 12.21 repeatedly. □ 


Now before ending this section and conversing solely of functions of exponential 
type, we see that functions like as well as are indefinitely summable. We 
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can find more complicated functions of order greater than one that are indefinitely 
summable if they behave nicely as the imaginary argument grows. 

3. The exponential space E and the indefinite sum 

As we’ve defined the indefinite sum at this point, we’ve only used functions 
bounded by an exponential as the imaginary part grows. We now consider bounding 
the function by an exponential as the real part of z grows as well. For this reason, we 
define an exponential space in which the indefinite sum operates. The exponential 
bounds will ensure we have a unique operator on the space and will show our second 
representation as well as our first representation are equivalent where they intersect 
in definition. 

Definition 3.1. Suppose f{z) is a holomorphic function on the right half 
plane z gC, 3?(z) > 0 and let | limjf( 2 )_>,o f{z)\ < oo. Assume there exists a,p,C G 
K+, 0 < q; < 7r/2, such that \f{z)\ < The space £ contains all 

such functions / that satisfy these conditions. 

£ contains polynomials, exponentials with a base that have a real part greater 
than zero, rationale with poles in the left half plane, and other special functions. 
Aside, on notation, we will write aj for the bound of / as the imaginary part grows 
and we will write pf for the bound of / as the real part grows (when distinction 
between functions is required). For our purposes we will be using the indefinite 
sum we defined in Theorem l3.ll Let / belong to £ and let cr G K'*', cr > 0, 


1 roo pt X pcT+ioo 

(3.1) d^dxdt 

2 ^ \^) Jo Jo Ja-ioo 

We know that if |/(z)| < that this representation converges and 

that I Yl,z /I < Therefore ^ is defined for any element of £. We also 

have the convenient formula, 

n 


which follows by the defining recursion of ^ and because / = 0 since lim^-^o ^ 


e^g{t) dt = 0. This limit is allowed to be taken because |e ^ e*g{t) dt\ < 


x—0 


CcX^ for any e > 0 as a; —?> 0 so that the integral expression for the differintegral 
converges at z = 0. We show now that ^ sends £ ^ £. 


X—0 


Theorem 3.1. Suppose f G £, then f G £. 

Proof. Let us say |/| < We know that / is bounded by 

an exponential and of type less than tt/2 on the imaginary line. Now observe that. 


E/= E /+E 

2 z-L5R(z)J j=l 


fo dt 
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which comes from the defining property of ^ and shows that, 


LSR(z)j 

IE/I< E ^ga|9(z)|+p|K(z)-LK(z)J+i| ^ M\\^{z)\ + 

z i=0 


And therefore since |[5ft(z)J + 1| < for some e > 0 this shows the 

result. 


□ 


With this we are now able to talk freely of the indefinite sum acting on a 
space. We can also give a more compact form for the operator which follows by 
Ramanujan’s master theorem. 


Corollary 3.1. Let f G E and define the entire function ^{x) = /O)) 

Then for z G <C, > 0, 


E/ = 


1 


n—0 i — 1 ^ ' 


ni-z) 


With these results we are prepared to analyze the indefinite sum in this space 
more intricately. 


4. On Indefinite Sum Convolution 

Beginning this section we will prove a lemma that appears elsewhere but we 
show it using techniques from fractional calculus. This result will be applied re¬ 
peatedly on various transforms involving the indefinite sum which will give us a 
powerful way of proving properties of the indefinite sum over the complex plane by 
looking at its values on the positive integers. 

Lemma 4.1. Let f G £ and for n gN assume /(n) = 0, then / = 0. 

Proof. Take 0 < ct < 1 and define g{x) = ^ ^(0/(1 - d^. By 

a simple exercise in contour integration we have g{x) = /(^ + 

implies that / = 0. □ 

We see this result follows quickly, however its power is abundant. Firstly, it 
implies if /i(n) = / 2 (n) for n G N then fi{z) = f 2 {z) for all 0 G We will 

not waste time in putting it to use. With this we define a new convolution that is 
associative, commutative, and spreads across addition. To be clear, we will make 
the following statement on notation, Xz / = Xz/('®)Vs, where the additional 
terms declare the variable we are performing the sum across. This notation allows 
us to nest more variables which is required for what comes next. 

Definition 4.1. Suppose f{z — l),g{z — 1) G £ and a/ -I- 2ag < 7r/2, then 
f ^ 9 = - 1 ) 5(2 + 1 - s) Vs 

z+l 
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We note this is definitively a convergent expression if we fix 2 ; and take the 
indefinite sum from Theorem O The motivation for introducing this convolution 
is because (/ x g){n) = J2'j=o ~ j) which is quite a familiar expression. 

Furthermore we know that the space £ is almost completely determined by how its 
functions behave on the naturals so this is bound to be valuable in some sense. We 
show some properties of this convolution acting in the space £. 

Theorem 4.1. Suppose f,g,hG £ then: 

(1) If af + ‘IoLg < 7r/2, then if f x g = h we have h€ £. 

(2) If af + 2ag < tt 12 and 2af + ag < 12, then f x g = g x. f. 

(3) If af + 2{ag + 2ah) < 7r/2 and then f x (g x h) = (f x g) x h. 

(4) If a f -\-2ag < 7r/2 and af -\-2ah < 7r/2 then f x (g + h) = {fxg) + {fxh). 

(5) If af + 2ag < tx 12 then f x g = 0 implies that / = 0 or g = 0. 

(6) Let {z + l)s_i = for s G C 5 f(^)>o and si,S 2 G Ck( 2 )>o we have 

r(J7) + l)si-i X r(^(^ + l)s 2 -i = r(si+s 2 ) l)si+s 2 -i- 

Prooe. We prove each result one by one. 

(1) First take f x g = Ez-LSR(^)J /(« “ + 1 - s) Vs + /(^ “ 

[5R(2)J — 1 + j)g([5R(2:)J + 1 — j). Therefore since |/(s — I)g{z + 1 — s)| < 

^ga^|9(s)|+a,(|a(s)| + |9(z)|) + (py:+p,)|Sf(s)| ggg | j _ 

s) Vs I < Where we treated 

as a constant independent of the variable s we are summing across pulled 

through the integral. 

(2) For the second, observe f x g and g x f are defined and if n G N, that 
(/ X g){n) = {g X f){n), so the result follows by Lemma STJ 

(3) For the third, by Lemma [4.II the result follows. 

(4) For the fourth, we observe the linearity of the indefinite sum and multi- 
pication. 

(5) For the fifth, we use Lemma Td.II again. As in X]j=o /(j)5(’^ ~ j) = 0 for 
all n implies f(n) = 0 or g{n) = 0 for all n. 

(6) For the sixth case, observe that /x 1 = X^z+i /(s—1) Vs so that 

l)„_i = 1 X 1 X 1 X ...{ntimes)... x 1. By associativity the result follows 
for integers so that + l)„-i x + l)m-i = + 

l)n+m-i- Now apply Lemma 14.11 on si and S 2 since G f in 

si,S 2 when 5R(0), 5R(si), 5R(s2) > 0. The result follows from this. 

This shows all the cases. □ 

With this we are prepared to show a new representation of complex iterations 
of the indefinite sum. To this end we define a new operator called the differsum. 
This term is novel to us and is used to express the operators similarity to the 
differintegral. 

Definition 4.2. Let (j){z) be holomorphic in the open set fl. A differsum Vz 
of (j) is holomorphic in s on open £1' and holomorphic in 2 : for z G fl, satisfying: 
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(1) For n S N we have = Vz Vz ■ ■ • {n times) ... Vz </> and Vz = 

(2) For so,si,so + si S n', Vz° Vz' </' = 

The definitive properties of the differsum is what we would expect of an iterate 
of the backwards difference operator/indefinite sum. It appears as an object that 
would be very difficult to analyze. Thankfully though, the deep connection between 
recursion and our differintegral allows us to speak freely of the differsum. 

Theorem 4.2. Suppose f{z) G £ and q G C, ^{q) > 0, then \/~‘^f{z) = 

J2z /(a)(^ + 1 — s)?-! Vs- 

Proof. We make use of the fact \y~^if{z — 1) = p^(/(-2) x (z + l)g_i)(z), 
where this is well defined because f G £ and because {z + l)g-i G £ in z and as the 
imaginary part of z grows it is bounded by for any e > 0. The proof of this 

is by induction for q an integer. It is clearly true for q = n = 1. Assume for n. Since 
j^^(z+l)n-i X 1 = :^(z+l)n and since = Vz+ifi^-^) x 1 = 

(/( 2 :)x-(;r^(z+l)n-i)xl = /(z)x(^;^((z+l)„_ixl)) = ffi/(z)x(z+l)„. Now 
since, for qi,q 2 G C, we have j^(z + l),,_i x p^(z + l)q,_i = r(g,+g,) (^) 9 i+g 2 -i 

we know that V^+i VJ”^^ ~ 1) = “ !)• Therefore removing the 

transfers up in z we have for 5R(q) > 0, 

vrv(s) = ^ 1 Vs 

□ 

To close we make use of our expression for the differsum. We take the transform 
on log(l+s). This will produce a generalization of the Gamma function. This comes 
about in the following sense, X]zl°§(l + a) Vs = log(r(l + z)), which follows by 
Lemma l4.ll Exponentiate and we have a well defined expression for F. Furthermore 
the iterates of the indefinite sum once exponentiated will produce functions ro(z) = 
z,Fi(z) = F(z),F 2 (z),F 3 (z), ... such that r„(z)r„+i(z) = r„+i(z + 1). These 
functions are normalized at one so that F„(l) = 1. 

We wish to be more brazen than this however and to develop the complex 
iterates of this recursion. Namely for z,q G C we have if 5R(g) > 0 and 5R(z) > 1 
that the function 

Fq(z) = iog(i+s)(z-s),_i Vn 

satisfies the recursion Fq(z)Fq+i(z) = rq+i(z + 1). This function expands the 
F function considerably and allows for a broad generalization of the Barnes G- 
function. In spirit, r 2 satisfies the same recursion as the Barnes G-function-we are 
unsure if the two agree however. 

5. Final Remarks 

We close hoping the reader has seen the connection between fractional calculus 
and recursion. The indefinite sum behaves quite neatly and we have only scratched 









ON THE INDEFINITE SUM IN FRACTIONAL CALCULUS 


13 


the surface of its structure and uses. We state that this paper is sister to [6]-on the 
differintegral and complex iterations. The two problems appear as an intrigue in it¬ 
eration and show the power of fractional calculus and its familiarity and connection 
with iteration in general. For a more clear application of this, observe [^. We are 
satisfied with the brevity of this paper as it is meant to show a pure mathematical 
curiousity which follows through fractional calculus. 

Although only briefly mentioned, our generalization of the Gamma function 
offers many questions. Can this functions be analytically continued to a mero- 
morphic function for q G C and z G C? What is its Weierstrass factorization in 
accordance to its multiplicative inverses zeroes? Explicitly, what is the Weierstrass 
factorization of p \ , in z dependent on q? 

i q \^Z) 

We are aware of more advanced techniques that apply to more complicated 
recursions. Many linear algebraic recursions in Tf{z) = f{z + 1) can be solved 
by these methods. Notably we solved for complex values of the recursion ■ 

We can do this for more complicated recursions by creating a similar isomorphisms 
between spaces that the Twisted Mellin transform induces. We can also perform 
fractional iterations of the recursion in a manner similar to how we produced the 
differsum. These results are very exotic and further capitulate the importance of 
the differintegral. 
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